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Abstract 

We find some errors in previous calculation of leading log QCD corrections to b —* 57 
decay, which include corrections from mt op to Mw in addition to corrections from Mw to 
mfc . The inclusive decay rate is found to be enhanced more than previous calculations. 
At nit = 170GeV, the running from mt op to Mw results in 13% enhancement, and for 
m t = 250GeV, 16% is found. 
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It is well known that the process b — > 57 is extremely sensitive to new physics beyond the 
Standard Model. In 1993, the CLEO collaboration placed an upper limit on the inclusive b — > 37 
decay B{b — > sj) < 5.4 x 10~ 4 at 95% C.L.[[L|. This has inspired a large number of papers[0. 
It has been argued that this experiment provides more information about restrictions on the 
Standard Model, 2-Higgs doublet model, Supersymmetry, Technicolor and etc. Recently CLEO 
has measured the inclusive branching ratio to be|| 

Br{b -> s 7 ) = (2.32 ± 0.51 ± 0.29 ± 0.32) x 10" 4 . (1) 

Corresponding to 95% confidence level, the range is 1 x 10~ 4 < Br(b — > 57) < 4 x 10^ 4 . More 
stringent constraints are obtained by experiments, so more accurate theoretical calculation of 
this decay rate is needed. 

The radiative b quark decay has already been calculated in many papers^, [5|, |7j. It is found 
to be strongly QCD-enhanced. In other words, the strong interaction plays an important role 
in this decay. However, there are still some uncertainties in these papers. Most papers |4], Q do 
not include the QCD running from m top to M^r. Since the top quark is found to be heavier than 
W boson( m top = 174 ± lOljJ GeV ), it needs a detailed calculation of this effect. Ref.0 
does include this running, however there are some errors in calculation of anomalous dimensions 
which may lead to some changes in final result. 

In our present paper, by using effective field theory formalism |J, we recalculate the b — *■ sj 
decay in Minimal Standard Model. We first integrate out the top quark, generating an effective 
five-quark theory. By using the renormalization group equation, we run the effective field theory 
down to the W-scale, including QCD corrections from m top to Myy, and correct some errors in 
ref.jJTJ. Then the weak bosons are removed. Untruncated anomalous dimensions of QCD running 
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from M\y to m b are used. Finally we calculate the rate of radiative b decay. 
The effective Hamiltonian is written as 



Heff = 2V2G F V tb V* Ci{n)Oiiji). (2) 



The operators Oi make a complete basis of dimension-6 operators: 



0\ R = -j^m b s L D 2 b R , 



2 LR = ^ 2 ^m b s L a^X% R G%, 

3 LR = ^ e/2 ^m b s L a^b R F, u , 
lo7r 

Qlr = /J, e glm b (f) + (f)_s L b R , 
P X l A = -^K a D»D»D°b L , 
Pi = ^^SLl^dTF^ 

p3 _ e/2 e Qb Tp „Mf) u h 
ID7T 

= i>V> e gi<t>+<l>-SL ]Z>b L , 
R\ = t^gl(D^ + )^s Lla b L , 

R 3 L = t^glMD a <P-)sLlab L . (3) 

The coefficients Cj(/i = m t ) are calculated from matching diagrams, and agree with ref.j7|. 
After evaluating the loop diagrams, we find that the weak mixing of operators agrees with 
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ref.0. While the QCD anomalous dimensions for each of the operators in our basis are 
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Comparing with ref.@, there are some differences in the anomalous dimension matrix, which 
may lie in omiting a symmetric factor of 1/2 in ref. in calculating Feynman diagram like Fig.l. 
And some changes may due to miscalculation. 

After these changes, the whole matrix can be easily diagonalized, and gives all real eigenvalues 
while that in ref. can not. Inserting anomalous dimension (|])([5D to the renormalization group 
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equation satisfied by the coefficient functions Cj(/x) 



d 



(6) 



we can have the coefficients of operators at fi = My/. And some of these operators change a lot 
from ref. 0j due to our improvements. 

In order to continue running the basis operator coefficients down to lower scales, one must 
integrate out the weak gauge bosons and would-be Goldstone bosons at /J, — Mw scale. This 
leads to the well-known six four-quark operators^, ||. The other remaining two-quark operators 
can be reduced by using equations of motion(EOM) to the gluon and photon magnetic moment 
operators 2 LR and 0\ R . 



To be comparable with previous results, we rewrite our operators 0\ R , 0\ R as O7, 0$ like 



ref-!, 



{e/16n 2 )m b s L a^b R F, 



8 = {g/l^ 2 )m b s L a^T a b R G 



(7) 



For completeness, we give the explicit expressions of the coefficient of operator Og and O7 at 



/! = M w , 
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They are expressed by coefficients of operators at ji = m t and QCD coupling a s . So it is 
convenient to utilize these formula. 

The obvious differences from QCD correction to C-j(Mw) and C%(Mw) can easily be seen 
from Fig. 2. In comparison to ref.0, the enhancement of coefficient of operator O7 is almost the 
same size, but the values for 0§ are quite different. Here the effect to Og is an enhancement 
rather than a suppression as in ref.jTj. These changes come from the corrections of anomalous 
dimensions described earlier. Since Cj(Mw) and C$(Mw) are both the input of the following 
QCD running from Mw to m b , it is expected to change the final result. 

The running of the coefficients of operators from \i = Mw to // = m b was well described in 
ref.||. After this running we have the coefficients of operators at fi = m b scale. Here we use 
M w = 80.22GeV, m b = 4.9GeV. Both C7(m b ) and C 8 (m b ) are enhanced in comparison to values 
obtained by Misiak where the QCD corrections from Mw to mt are neglected ||. 

The leading order b —> sj matrix element of H e ff is given by the sum of operators O5, 0§ 
and O7, this disagrees with ref.0, but agrees with Misiak||]. The sought amplitude will be 
proportional to the squared modulus of 

C e 7 ff (m b ) = C 7 (m b ) + Q d [C 5 (m b ) + 3C 6 (m b )] (10) 

instead of \C 7 (M b )\ 2 itself. 

Following ref. |5|, applying eqs.(pH|), one finds 

BR(B^X sl ) ^ 6a QED ^eff^ (l- 2 < m ^firn c / mb )\\ (11) 



BR(B — > X c eu) Tig{m c /m b ) \ 37T 

where g(m c /m b ) ~ 0.45 and f(m c /m b ) ~ 2.4 corresponding to the phase space factor and the 
one-loop QCD correction to the semileptonic decay, respectively |5|]. The electromagnetic fine 
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structure constant evaluated at the b quark scale takes value as aQED{mb) = 1/132.7. The 
results are summarized in Fig.3 as functions of the top quark mass. The QCD-uncorrected 
values are also shown. In this figure, one can easily see that, at m t = 170GeV, it results in 13% 
enhancement from Misiak's result 0, and for m t = 250GeV, 16% is found. 

As a conclusion, we have given the full leading log QCD corrections(include QCD runnings 
from rritop to Mw), with whole anomalous dimension matrix untruncated. Comparison to the 
previous calculation j7|, two points are improved: 

(1) Correct errors of anomalous dimensions in ref. [^]. 

(2) Use untruncated anomalous dimensions in QCD running from M w to m b instead of 
truncated ones. 

In fact, point(l) makes an enhancement while point(2) leads to a suppression. The total 
result does not change a lot, e.g. a suppression of 3% at m top = 170GeV comparing ref. 0. 

The whole Q CD-enhancement of the BR(B — > X^) makes a factor of 3.9 at m t = 170GeV, 
when Aqcd = 175MeV. Using the experimental branching ratio BR(B — > X c eV) = 10.8%, one 
can find that BR(B — > Xg'-f) ~ 4 x 10~ 4 at m t = 174GeV. It just reaches the upper limit of 
present experiment of CLEO. That shows there is very little space for new physics. 
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Figure Captions 

Fig.l One of the Feynman diagram in calculating Anomalous dimensions, with the heavy dot 
denoting high dimension operator. 
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Fig. 2 The photon and gluon magnetic moment operator's coefficient Cj(Mw) (upper) and 
C 8 (M W ) (lower) for different top quark mass. The ones with and without QCD corrections are 
indicated by solid and dashed lines respectively. (Aq C d = 300MeV is used) 

Fig. 3 BR(B — > x s j) normalized to BR(_B — > x c eV), as function of top quark mass. The upper 
solid lines indicated our results for a full QCD correction. Dashed lines correspond to Misiak's 
results without QCD running from m top to M w . 
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